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ABSTRACT 

The UV spectra of nova-like variables are dominated by emission from the accretion disk, modified by 
scattering in a wind emanating from the disk. Here, we model the spectra of RW Tri and UX UMa, the only 
two eclipsing nova-like which have been observed with the Hubble Space Telescope in the far-ultraviolet, in an 
attempt to constrain the geometry and the ionization structure of their winds. Using our Monte Carlo radiative 
transfer code, we computed spectra for simply parameterized axisymmetric biconical outflow models and were 
able to find plausible models for both systems. These reproduce the primary UV resonance lines - N V, Si IV, 
and C IV - in the observed spectra in and out of eclipse. The distribution of these ions in the wind models is 
similar in both cases as is the extent of the primary scattering regions in which these lines are formed. The 
inferred mass-loss rates are 6% - 8% of the mass accretion rates for the systems. We discuss the implication of 
our point models for our understanding of accretion disk winds in cataclysmic variables. 
Subject headings: accretion, accretion disks — binaries: close — novae, cataclysmic variables — radiative 
transfer — stars: individual (RW Tri, UX UMa) — stars: winds, outflow 



1. INTRODUCTION 

Mass outflows are associated with accretion via a disk onto 
compact objects in a wide range of astrophysical systems, 
ranging from proto-stellar objects to close binaries and ac- 
tive galactic nuclei. The presence of mass outflow in all these 
systems, which cover a wide range of scales and system pa- 
rameters, suggests an important connection between the pro- 
cesses of accretion and the outflow. Cataclysmic variables 
(CVs) are the closest astrophysical examples of such systems, 
and hence are an important test case for our understanding 
of outflows from disks. CVs are close binary systems con- 
sisting of a white dwarf (WD) and a low-mass late-type sec- 
ondary star The secondary overflows its Roche lobe, lead- 
ing to mass transfer to the WD. Due to angular momentum 
conservation, the transferred material builds up an accretion 
disk around the WD (unless the magnetic field of the WD is 
large enough to disrupt the flow). Disk emission usually dom- 
inates other sources of light at visible wavelengths. Many 
CVs undergo outbursts of 3 - 5 mag lasting from days to 
weeks on timescales of weeks to years. The outbursts are due 
to a thermal instability in the disk which transforms the disk 
from a relatively cold, largely neutral gas to a hot, ionized 
plasma, increasing the accretion rate in the inner dis k from 
about 1 X IQ-" Mp, yr-^ to 1 0"^ - 10"^ Mq yr"! llRshil 
11979: Mineshig e & Osaki ^ 1983). The frequency of outbursts 
depends on the longer term mass-transfer rate from the sec- 
ondary, which for systems that underg o outbursts mus t be be- 
tween the quiescent and outburst rate ('Canni zzolll993h . 

Far-ultraviolet (FUV) spectra of CVs obtained with the In- 
ternational Ultraviolet Explorer (lUE) were the first to show 
that CVs drive winds when the disk is in the high mass- 
transfer sta te (Ki-autter et al. 198U iGreenstein & OkdlT982l 
ICordova & Masoa.1982) . Winds are not observed when CVs 

' Space Telescope Science Institute, 3700 San Martin Drive. Baltimore, 
MD 21218, USA 

-Max Planck Institute for Astrophysics, Karl-Schwarzschild-Str. 1, 
85748 Garching, Gemiany 

^ Department of Physics, University of Southampton, Southampton 
S017 IBJUK 



are in the quiescent, low mass-transfer state. The signa- 
ture of the wind is imprinted on the spectra in the form of 
P Cygni-like and/or blueshifted absorption profiles of res- 
onance lines, such as N V, Si IV, and C IV. The blue 
edges of the absorption profiles imply velocities exceeding 
3000 kms~^ in s ome systems (see, for a recent example, 
[Hartley et alj l2002h . The lUE observations were used to 
show that the emission features are generaUy weak or ab- 
sent in systems viewed at low inclin ation and progr essively 
more prominent at higher inclination (lLaDouslll991h . Time- 
resolved lUE observations also demonstrated that line shapes 
in eclipsing and nearly eclipsing systems vary systematically 
through secondary conjunction, showing not only that the 
emission lines had to arise from a region at least compara- 
ble in size to that of the secondary star but also that the out- 
flowing material was likely rotating (Cordova & Mason 1985]; 
iDrew & Verbunt 1988). While initial modeli ng of CV winds 
favor e d spherical geometry for s implicity ( Drew & VerbuntI 
119851: iMauche & Raymondlll987l) . the observed characteris- 
tics of the resonance lines in CVs have over time led to a con- 
sensus picture of a biconical mass outflow originating from 

th e inner disk. 

IShlosman & Vitellol (119931) . hereafter SV93, were the first 
to develop a radiative transfer code based on the Sobolev ap- 
proximation that was capable of determining the ionization 
structure of a biconical outflow emerging from the disk and 
of modeling C IV line profiles. They then used this to at- 
tempt to constrain the wind geometry of non-magnetic CVs. 
IVitello & Shlosman (1993) found that they could reproduce 
the observed profiles of C IV in three systems — RW Sex, RW 
Tri, and V Sge — with plausible geometries and mass-loss 
rates; they also found that they could produce similar C IV 
profiles with quite different accretion rates by changing other 
parameter s in their kinematic description of the wind. Sub- 
sequently.lKni gge. Wood & Drew I (Il995h . hereafter KWD95, 
developed another Monte Carlo radiative transfer code which 
solved the radiative transfer exactly using a slightly different 
kinematic prescription for a biconical wind. Their code did 
not calculate the ionization balance in the outflow and instead 



Noebauer et al. 



assumed constant ionizati on fractions throughout. As had 
IShlosman & Vitellol (119931) . they were able to reproduce qual- 
itatively many of the basic characteristics of C IV line profiles 
in disk-dominated CVs. Knigge & Drew ( 1997) applied their 
code to high-resolution time-resolved Hubble Space Tele- 
scope (HST) spectra of C IV of the eclipsing nova-like vari- 
able UX UMa, concluding that one not only needed a fast 
wind with a wide opening angle but a relatively dense, slow- 
moving transition region between the photosphere of the disk 
and the wind. 

Despite the clear observational indications for the pres- 
ence of mass outflow in CVs, much is still unknown about 
the physics of these winds, in particular the driving mecha- 
nisms. By analogy with O-star winds, most have assumed 
that the winds of CVs are radiatively driven. Hydrodynamical 
simulations seem to broadly support this picture (Proga et aT] 
12002), even though questions regarding whether CV disks are 
luminous enough to d rive the required mass outflow remain 
jDrew & Progall20^0l) . Alternatives, magnetically driven or 
magnetically assisted flows from the disk or boundary layer 
region exist, base d on suggestions tha t ca n be traced to , 
among others. iBlandford & Paynl (1 1 982h and lPringld (Il989h . 
and are widely discussed for other disk systems, ranging 
fro m young stellar ob j ects to act ive galactic nuclei (see, e. 
g.. lKoniglet"an 120101: ICampbelll 12010.) . An important ob- 
stacle to a confident identification of the driving mechanism 
is that even basic wind parameters, such as the mass-loss 
rate, the acceleration, collimation, and ionization structure, 
remain poorly constrained. Thus obtaining reliable informa- 
tion on the wind structure, and especially about the ionization 
in the wind, should lead to a much better understanding of the 
physics underlying CV winds. In the long term, this might 
even allow us to gain insight into other types of accreting sys- 
tems, as mass outflows seem to be a general feature of such 
systems. 

Nova-like variables are disk-dominated CVs, which, unlike 
most other CVs, are nearly always observed in the "outburst" 
state. This is presumably because they have a higher mean 
mass-transfer rate than othe r types of CV (e.g., dwarf novae; 
ISchreiber & Gansickel2002h . Nova-like variables have orbital 
periods just above the period gap in non-magnetic CVs, where 
magnetic braking is thought to produce the largest mass- 
transf er rates in the long-term evolution of CVs JHowell et alj 
1200 Ih . Because they stay in the high mass-transfer state for 
long periods of time, they may be the closest examples to 
steady-state accretion disks among the CVs. RW Tri and UX 
UMa are the two brightest eclipsing examples of this class. 
Eclipsing systems are of particular interest in the study of 
CVs, and, in particular, the winds of CVs, since time-resolved 
observations of the eclipse provide spatial information about 
the disk and wind properties that is unavailable for systems 
viewed at lower inclinations. 

Here, we describe our efforts to model the HST spectra of 
RW Tri and UX UMa in an attempt to place constraints on 
the wind geometry and the ionization structure. These are 
the only eclipsing nova-like CVs for which medium spec- 
tral resolution, time-resolved UV spectra exist. The obser- 
vations of RW Tri have been briefly discussed by Mason et alj 
( 11997) ; the observations of UX UMa have been discusse d 
more extensively by^^ ason et al. ( 1995), Ma son et alJ(ll997h. 
iKnigge & Drew! (f!997.). .Baptista et al., a998l) . iFroning et alj 
( I2()03h. and iLinnellet all ( 120081) . \\^h the partial exception 
of iKnigge & Drew! (119971) . the earlier studies were primar- 
ily focused on modeling the accretion disk and WD in UX 
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Figure 1. Light curves of the two RW Tri observations that were used to con- 
struct average pre-eclipse and mid-eclipse spectra. The phase of each mea- 
surement was derived with the revised ephemeris of RW Tri of Groot et al] 
I2004i) . The measurements in the upper panel all cover the shorter wavelength 
band from 1 148 A to 1425 A, the observations in the lower panel were taken 
in the band from 1376 A to 1663 A. The phase intervals that were used for 
the averaging process are shaded in gray. 

UMa. However, no one has attempted to constrain the prop- 
erties of the wind in these systems by modeling all of the 
prominent wind lines simultaneously. The remainder of this 
paper is organized as follows. In Section |2l we describe the 
time-resolved observations of RW Tri and UX UMa with the 
Goddard High Resolution Spectrograph (GHRS) and the se- 
lection of the spectra we attempt to model. In Section [3] 
we describe briefly the radiative transfer code developed by 
iLong & Kniggd ( 120021) to produce synthetic spectra from pa- 
rameterized models of biconical outflows, which we use to 
generate synthetic spectra. Then, in Section HI we describe 
the strategy we used to identify models for the systems, and 
in Section |5] we discuss the degree to which we were able to 
reproduce the observe spectra and the implications of the re- 
sults on the wind geometry of RW Tri and UX UMa. Finally, 
in Section|6l we summarize what we have learned through our 
attempt to model the pre-eclipse and eclipse spectra of RW Tri 
and UX UMa. 

2. OBSERVATIONS 

UV spectra of the nova-like RW Tri and UX UMa were ob- 
tained with the GHRS instrument on HST in late 1996 and 
early 1997 (as part of programs 6494 and 6024) and are de- 
scribed by Mason et al. ( 1997). Both systems were observed 
in order to explore time variations in the spectra as a function 
of the orbital phase. All of the observations discussed here 
were performed with the G140L grating, providing a wave- 
length resolution of 0.57 A, in the rapid read-out mode, in 
which spectra were recorded at 5 s intervals. Each observa- 
tion covers one of two overlapping wavelength intervals, each 
with a bandwidth of about 287 A so that spectra from sepa- 
rate observations must be combined to yield the full spectral 
range from 1 150 A to 1660 A. The observations of UX UMa 
and RW Tri each comprised about nine HST orbits. An ob- 
servation log is presented in Table [T] To carry out this study, 
we retrieved and used the final reprocessed data from the HST 
archive. 

To facilitate comparison of the synthetic spectra to the ac- 
tual spectra of RW Tri and UX UMa, we averaged the ob- 
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Table 1 

Observation Log for HST Spectra of RW Tri and UX UMa. 

ObsID Obs. Start (MJD) Spectral range (A) Pliase Coverage 
UXUMa 



z3fy0304t 

z3fy0306t 

z3fyal02t * 

z3fya302t 

z3fybl02t 

z3fybl04t 

z3fybl06t 

z3fybl08t 

z3fyblOat " 



50398.2014 
50398.2132 
50395.1085 
50398.2554 
50395.1751 
50395.1998 
50395.2422 
50395.2526 
50395.3091 



1148-1435 
1376-1663 
1148-1435 
1148-1435 
1376-1663 
1148-1435 
1148-1435 
1376-1663 
1376-1663 



0.654 ■ 
0.714- 
-0.073 ■ 
-0.071 ■ 
0.027 ■ 
0.391 ■ 
0.607 ■ 
0.660 ■ 
0.053 ■ 



0.701 
0.762 
0.033 
0.093 
0.377 
0.431 
0.647 
0.771 
0.111 



RWTri 



z3kxal02t 

z3kxal04t 

z3kxal06t * 

z3kxa202t 

z3kxa204t 

z3kxa302t 

z3kxa304t 

z3kxa306t * 

z3kxa402t 

z3kxa404t 

z3kxbl02t 

z3kxb302t 



50471.3006 
50471.3446 
50471.4116 
50470.6978 
50470.7418 
50474.3154 
50474.3593 
50474.4263 
50473.7126 
50473.7566 
50471.4786 
50474.4933 



1148-1435 
1148-1435 
1148-1435 
1148-1435 
1148-1435 
1148-1435 
1376-1663 
1376-1663 
1376-1663 
1376-1663 
1376-1663 
1376-1663 



0.395 ■ 
0.585 ■ 
-0.126. 
0.796 ■ 
-0.015 ■ 
0.395 ■ 
0.585 ■ 
-0.126. 
0.796 ■ 
-0.015 ■ 
0.162- 
0.163- 



0.433 
0.718 
0.007 
0.893 
0.118 
0.433 
0.718 
0.007 
0.833 
0.118 
0.296 
0.296 



Note. Average spectra are constructed from observations labeled with an 
asterisk (*). 



served spectra over specific phase intervals. There were cer- 
tain difficuhies associated with this. At any time observa- 
tions only covered parts of the spectral range of interest and 
the phase coverage was also not ideal for either UX UMa or 
RW Tri, as indicated in Table 1 . For RW Tri, two observa- 
tions sampled pre-eclipse and mid-eclipse orbital phase. The 
other observations covered only parts of the mid-eclipse and 
the post-eclipse state. In both cases, one observation existed 
for each wavelength band. For UX UMa, there was only one 
observation covering the eclipse in the 1375 - 1665 A band, 
but two covering the shorter wavelength band from 1148 A 
to 1435 A. To construct representative spectra for the mid- 
eclipse and the pre-eclipse phase, both RW Tri observations 
containing information of these two states were taken. For 
UX UMa, the two observations in the 1148 - 1435 A band had 
different overall pre-eclipse flux levels but showed no striking 
differences in spectral shape. Therefore, only the observation 
which was closest to the overall flux level of the only measure- 
ment set in the longer wavelength band was used. The spectra 
used in our study were obtained by averaging over represen- 
tative phase intervals. Choosing the width and the position of 
these intervals was, at least in the pre-eclipse case, somewhat 
arbitrary, but tests showed that the overall flux level and the 
spectral features in the average spectra were not very sensitive 
to these choices. In the end, we used the highlighted regions 
in Figures [T] and |2] for the averaging process. The resulting 
spectra were finally merged to cover the entire wavelength 
band from 1 148 A to 1665 A by averaging over the overlap- 
ping region. Figures |3] and H] show the resulting pre-eclipse 
and eclipse spe ctra for UX UMa and RW Tri, respectively. 
According to .Rutten et alj ( 11992b . RW Tri is thought to lie 
along a line of sight with E(B-V) of 0.1, and so to facilitate 
comparisons with models the spectra for RW Tri were cor- 
rected for reddening. Reddening along the line of sight to UX 
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Figure 2. Light curves of UX UMa. The phase of e ach measurement wa s 
determined with the ephemeris of UX UMa derived by'Baptista et al.'('199y). 
The upper panel shows the observation in the short wavelength band. In the 
lower panel, the observation covering the long wavelength band is displayed. 
The average spectra were constructed from the shaded regions. 
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Figure 3. Average RW Tri HST pre-ecHpse and mid-eclipse spectra obtained 
by averaging over the regions highlighted in Figure [T] In the pre-eclipse 
panel, the result of fitting the continuum to an accretion disk radiating as an 
ensemble of stellar atmospheres is also shown. The main emission lines are 
shown in the lower panel. 



UMa is believed smaU (E(B-V}<0.04. IVerbunti[T987h . so no 
correction was made for UX UMa. 

The pre-eclipse spectra of both RW Tri and UX UMa show 
rather flat continua with a number of superimposed resonance 
lines. During eclipse the continuum is suppressed and the 
spectrum is dominated by the emission lines that, because the 
disk is eclipsed, must arise from an extended region. 

The three major resonance lines present in both systems 
are the doublets N V AA1239, 1243, Si IV AA1394, 1403, 
and C IV AA1548, 1551. Apart from these there are addi- 
tional, weaker lines of low-ionization state metals such as C II 
A1333, C III A1176, and Si III AA1295 - 1303 and lines 
of both hydrogen (H I Lyman a A1215) and helium (He II 
A1640) are visible. In the pre-eclipse state of both systems, 
the low ionization state lines, especially those of C II A 1333 
and Si III AA1295 — 1303 are seen as deep absorption dips. 
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Figure 4. Average UX UMa//Sr pre-eclipse and mid-eclipse spectra. Again 
emission lines are identified in the lower panel and the disk continuum fit is 
shown in the upper panel. 

suggesting that the main concentration of these ions lies in the 
line of sight. The UX UMa pre-eclipse spectra generally show 
many absorption features, including in the main resonance 
lines: for example, two prominent narrow dips are present 
in the Si IV line and one in the C IV line. 

3. RADIATIVE TRANSFER CODE AND WIND MODELS 

In order to place limits on the geometry of the winds in 
RW Tri and UX UMa using the spectra, we need to model 
the effect of the wind on the spectra. Our working hy- 
pothesis is that the winds arise in a biconical outflow. To 
model the spectra we parameterize the outflow using one 
of two kinematic descriptions for the wind and we use a 
Monte Carlo radiative transfer code (PYTHON, developed 
by Long & Kni gge 2002, with improvements described by 
ISim et all 120051) to calculate the ionization structure of the 
wind and to create synthetic spectra. The chemical compo- 
sition of the wind material is fixed to th e solar abundances 
determined bv ' Anders & Grevessd (I1989I) . Here, we give a 
short summary of the important features of this program and 
the two classes of parameterize d wind mode ls which we use 
in this study (for full details, see lLong & Knigge 2002) . 

3.1. Monte Carlo approach 

PYTHON operates by performing a sequence of Monte 
Carlo simulations in which the quanta are packets of radiative 
energy. These packets are launched in accordance with the 
specified radiation sources (see below) and propagate through 
a parameterized model for the CV disk wind. As the pack- 
ets propagate, they can interact with the wind material via 
electron scattering, resonance line scattering, bound-free ab- 
sorption, or free-free absorption. Their trajectories are finally 
terminated when they either escape, lose all their energy to ab- 
sorption in the wind, or are lost from the simulation by strik- 
ing the WD, disk surface, or companion star. 

Since the ionization and heating of the wind are controlled 
by the radiation, the ionization and thermal state must be com- 
puted self-consistently with the local radiation field proper- 
ties. This is achieved via an iteration sequence. In each 
iteration, the wind properties are held fixed while a Monte 
Carlo simulation of the radiation field is performed. During 
the Monte Carlo simulation, estimators for the local radiation 



field properties (radiation temperature Tr and dilution factor 
W) are recorded for each grid cell in the wind and used to 
make an improved estimate of the wind ionization state and 
kinetic temperature Tc for use in the next iteration. 

Once the wind properties reach convergence, a final se- 
quence of additional Monte Carlo simulations for the con- 
verged wind properties are used to obtain the observed spec- 
trum for the desired inclination angle. In these final simula- 
tions, a modified version of the technique described bv lWoodsl 
(119911) is used: in every interaction with the wind medium the 
packet contributes to the final spectrum according to the prob- 
ability of emission into the observer's line of sight after the 
interaction. 

3.1.1. Radiation sources 

The calculations presented in this work include radiation by 
the WD and its accretion disk. The disk is the primary source 
of radiation in non-magnetic CVs in the high mass-transfer 
state. It is treated as a geometrically thin but optically thick 
disk which follows the standard temperature distribution (e.g., 
I Wade 1984) 
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The disk radiation is composed from an ensemble of 300 
concentric annuli each of which radiates a spectrum at the 
local value of Toff. The inner and outer boundaries of each 
annulus are determined by the requirement that all annuli 
should emit an equal fraction of the total disk luminosity. 
For this study, we model the spectrum from each disk an- 
nulus as a black body during the iteration processes used 
to obtain the ionization state. In the spectrum calculation, 
the disk emission is described by synthetic stellar spectra 
created with TLUSTY/SYNSPEC dHubenv & Lanzl [l995.) . 
When synthetic spectra are used, the calculation of the ef- 
fective gravit y of t he disk follows the prescription given by 
Herter et al.i (119791) . To describe the WD radiation, a sin- 
gle black body/Hubeny model is used. The secondary star 
is treated as a dark absorber. 

The PYTHON code can also account for radiation from a 
boundary layer par ameterized by a lumino sity and tempera- 
ture as described bv lLong & Kniggj (120021) . We found, how- 
ever, that if we included a hot luminous boundary layer (such 
as introduced by Knigge & Drew 1997) the wind was typi- 
cally too highly ionized to obtain a good match with observed 
spectra. We therefore opted to exclude any such boundary 
layer from the calculations. This may imply that either the 
WD has significant rotation or that the boundary layer geom- 
etry is such that much of the wind is not exposed to its radia- 
tion. 

The energy emitted by the radiation sources is discretized 
into energy packets, each with its own direction of propaga- 
tion, frequency, and weight. The packets are initialized with 
weights corresponding to the fraction of the physical luminos- 
ity carried by the package and frequency sampled randomly 
from the luminosity spectra of the WD and accretion disk an- 
nuli. Initial directions of propagation are assigned assuming 
a linear limb-darkening law and initial positions are located 
on either the WD or disk surface, as appropriate. In addi- 
tion to the primary radiation sources, the code also accounts 
for free-free, free-bound, and bound-bound emission by the 
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wind, based on its local ionization state and kinetic tempera- 
ture. 

3.1.2. Wind conditions 

A key component of modeling the wind spectrum is a 
realistic treatment of the wind ionization state. In con- 
trast to some previous studies in which uniform ioniza- 
tion conditions throughout t he wind were assumed (e.g., 
iKnigge. Wood & Drew Il995h . python computes the ioniza- 
tion state as a function of position using the modified nebular 
approximation described by Mazzali & Lucy ( 1993]) 



^W[^ + W{1 



■ rj. N 1/2 



T'j + l'^c 



(2) 



Here, W is a dilution factor, ^ is the fraction of recombina- 
tions that go directly to the ground state, and Tc and Tr are the 
electron and radiation temperatures, respectively. The aster- 
isk denotes quantities to be evaluated in local thermodynamic 

equilibrium at 2^^; 

In contrast to IMazzah & Lucvl (Il993h . however, we do not 
assume a simple relationship between Tr and Te but instead 
calculate To independently from the assumption of thermal 
equilibrium (i.e., the heating and cooling rates are in equilib- 
rium at every point in the wind). As noted above, Tr, W, and 
Tc are not known a priori, but are determined by an iterative 
sequence of Monte Carlo experiments. In practice, Tr is the 
most important parameter for the ionization balance (given by 
Equation (|2]i) and is primarily determined by the temperature 
distribution in the disk. This in turn depends on the adopted 
system parameters (WD mass and radius, and the mass accre- 
tion rate M). 

3.2. Wind models 

For our simulations, we adopt two classes of simply param- 
eterized wind models, those developed by SV93 and KWD95. 
We carried out our studies using both wind models so that we 
could investigate any sensitivity of our results to the parame- 
terization. As demonstrated below, however, we find that our 
conclusions are relatively insensitive to this choice. Both of 
these models describe the wind as a stationary biconical out- 
flow emerging from both sides of the accretion disk. In each 
case, the velocity in the wind consists of a rotational compo- 
nent Vtj, and a poloidal component Wq. The rotational compo- 
nent is obtained assuming that the wind streamlines rise from 
a disk in Keplerian rotation and that specific angular momen- 
tum is conserved as the flow moves outward: 



v^{r)r 



"(0,0) 



{R)R, 



^(0,0) 



(i?) = 



GM. 



V\fD 



R 



(3) 



(4) 



Here, and in the following, r stands for the radial position of 
a point in the wind, while R denotes the radius at which the 
wind streamline passing through the point reaches the accre- 
tion disk. The poloidal velocity field in the wind causes mate- 
rial to spiral on streamline cones away from the disk. The de- 
tailed geometry of these cones depends on which of the wind 
prescriptions is used. This, together with the way of parame- 
terizing the poloidal velocity and calculating the wind density, 
is the main difference between the two classes of wind models 
we consider. 
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flmm fl,„„x 



Figure 5. Left: wind geometry as parameterized by KWD95. The black 
circle denotes the WD and the hoiizontal line the accretion disk. Note that 
the wind shape is controlled by only one collimation parameter, d, which 
is usually given in units of RwD- Right: wind geometry as parameterized 
by SV93. In both panels, the mass outflow region is shaded in gray. For 
illustrative purposes only one half of the wind structure is shown (the wind is 
symmetric under reflection in the xy-plane). 

3.2.1. KWD95 wind model 

KWD95 parameterize the wind as a biconical outflow from 
the entire accretion disk and characterize its geometry by only 
one parameter d. By definition, this parameter specifies the 
position of a point on the rotational axis where extrapolations 
of all streamline cones meet. Together with the disk extent 
[-RwDi ^disk]. d defines the opening angles of the wind. Fig- 
ure |5] illustrates the wind geometry of the KWD95 treatment. 

The wind region itself is characterized by a density and a 
velocity. The density in the wind is specified by the mass-loss 
rate per area in the disk. This quantity is assumed to follow 
the temperature distribution in the disk 



rniR) = M 



wind 



T^SiR) 



^dA'T^g{R') 



(5) 



The integral in the denominator ensures the normalization to 
the total mass-loss rate M^ind when integrated over the disk 
surface from which the wind emerges. The exponent a allows 
the mass-loss rate per unit area to be tied to the local disk 
effective temperature (a — 0.25) or the luminosity (a = 1.0) 
in the disk. The mass-loss rate per area can be used together 
with the poloidal velocity W(q,o) at the streamline base and the 
polar angle of the streamline 5 = arctan -j to calculate the 
density at the base of the streamline 



Po{R) 



m{R) 



V(qfi)iR)cOsS 

As the wind rises from the disk, the density decreases as 



p(r) = po{R) 



qcosSj Vq{r) 



(6) 



(7) 



Here q denotes the distance to the point d on the rotational 
axis. The poloidal velocity component Vq{T) depends on three 
parameters, the acceleration length Rs, acceleration exponent 
/3, and terminal velocity Wt»: 



Vqil) = Cs{R) + [Voo - Cs{R)] 1 



Rs 



Rs + l 



(8) 



The terminal velocity is parameterized as a multiple (/) of the 
local escape speed 
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fVc 



I 



2GAI^^-£, 



R 



(9) 



and the poloidal velocity is assumed to be equal to the local 
sound speed at the base of the wind streamlines. 



Cs(i?) = 10 



104 



km s 



(10) 



The quantity / in Equation (O represents the distance of the 
point r to the accretion disk, measured on the streamline cone. 
In total there are six parameters, d, a, /3, A/wind, /, and 
i?s which specify the CV wind in the KWD95 parameteri- 
zation. Although the true geometry is controlled only by d, 
the extent of the wind region in which the ion concentrations 
are sufficiently high to lead to strong spectral lines is strongly 
sensitive to the distribution of mass outflow from the disk, as 
determined by M^ind and a. In particular, high values of the 
exponent a will cause a very high density near the rotational 
center of the wind and a very dilute outflow at the outer edges 
of the wind while decreasing a moves more material to the 
outer edges. Thus, the effective coUimation is determined by 
a combination of d and a. 

3.2.2. SV93 wind model 

SV93 presented a wind model that is more flexible than the 
KWD95 treatment, but involves more free parameters. In the 
S V93 description, the wind does not necessarily emerge from 
all parts of the disk. Instead, the wind region is specified by a 
minimum and a maximum disk radius i?,nin and i?max- Sim- 
ilarly, the collimation of the wind is specified by a minimum 
and a maximum streamline angle 0nii„ and 0„iax. Figure |5] 
illustrates the wind geometry in the SV93 model. The distri- 
bution of streamline angles throughout the disk is determined 
by an additional parameter 7 defined via 



e{R) = 6, 



MR)' 



where 



xiR) = 



R — Ru 



-'T^max J^r\ 



(11) 



(12) 



As in the KWD95 description, R is the radius of the point at 
which the streamline passing through point r crosses the disk 
plane. 

The rotational velocity in the SV93 wind model is exactly 
the same as in the KWD95 description (see Equations (O and 
(|4|i). The poloidal velocity (vq) at a distance I measured on 
the streamline cone is specified in terms of three parameters: 
an acceleration length scale Rs, acceleration exponent a, and 
the terminal velocity Woo 



^q(0 = Wo + (Woo - Vq) 



(i/Rsr 



{l/Rs)'' + 1 



(13) 



We follow SV93 in adopting wq = 6 km s^^. As in the 
KWD95 prescription, the terminal velocity is specified in 
units of the local escape speed (see Equation (|9]l). 

Again, the wind density is determined by the disk mass- 
loss rate per unit area m. In the SV93 prescription, this is not 
specified in terms of the local disk temperature but is assumed 
to follow the relation 



rn{R) = Mwind 



R^cosejR) 
J dA'R'^ cos 9{R') 



(14) 



The normalization ensures that calculating the total mass-loss 
rate is given by A/^ind while A allows for variation in the mass 
loading of wind streamlines. The density in the wind is then 
given by 



p{r,z) 



RdR m{R) 
r dr Vz{r, z) 



(15) 



The R/r dR/dr factor results from the increasing area be- 
tween streamlines the further the wind is away from the disk 
which, in the special case of streamlines on a cone with angle 
9, simplifies to 



dr _ l(de/dR) 
dR cos 6 



(16) 



As in the KWD95 description, / is the distance to the disk on 
the streamline cone. 

In total, the SV93 prescription requires ten parameters to 
define the wind: 7, i?„ii„, i?niax, 6*^111, ^max. A, R^, a, /, and 
A^wind- The collimation of the wind is primarily determined 
by ^min and 0niax> Supplemented by the effects of A. 



4. FINDING FIDUCIAL MODELS: BASIC APPROACH 

Our primary objective is to constrain the physical properties 
of the wind. The empirical wind models described in the pre- 
vious section involve a rather large number of parameters, too 
many for it to be computationally practical to compute large 
grids of model spectra that explore the complete parameter 
space. Therefore, we have chosen to accept previous determi- 
nations of the system parameters, where possible, and focus 
our study on the major parameters of the wind itself. 

The distance to RW Tri has been determined to h& D — 
Ml'^ll pc from its paraUax dMcArthur et al.ll999h. Th i s is al- 
most identical to the distance obtained bv lRutten et al.l d 1 992b 
from the analysis of the optical eclipse light curve. Accord- 
ing to Rutten et all (11992) the WD mass for RW Tri is about 
0.7 A/0, which is typical of WDs in other systems and the 
mass ratio is fairly close to 1 . In their study they also estimate 
that the inclination for RW Tri is about 75°, although this is 
somewhat uncertain and v alues as low as 67° iKaitchuc k et al.l 
119831) and as high as 80° dMason et al."1997) have been sug- 
gested. From spectroscopic studies, Kaitchuck et al. (198^ 
also determined the K velocity of the WD. The amplitude of 
this velocity is roughly 200 km s^^. Table|2]lists the system 
parameters we adopted for RW Tri. 

The main source for the continuum flux in nova-like vari- 
ables is normally the accretion disk. For a given wind model, 
the ionization structure is primarily determined by the radia- 
tion temperature Tr, which, for a steady state, is controlled 
by the effective temperature {Tcfd of the disk. Assuming 
that the disk radiates as an ensemble of stellar atmospheres, 
as described by Equation ([T]i, then the mass accretion rate 
is a function of the observed flux, the distance, inclination, 
and the mass and radius of the WD. Thus, we can use the 
observed continuum flux and system parameters to constrain 
the mass accretion rate and avoid treating it as a completely 
free parameter in our studies. For this exercise, we restricted 
the wavelength range for the continuum fits to the regions 
between the Si IV, C IV, and He II lines which show a 
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Table 2 






System 


Parameters for RW Tri and UX UMa. 


System 


Parameter 


Value 


Reference 


RWTri 


AfwD 


0.7 Af0 


Ruttenetal. (1992) 




Afa 


0.6 Af0 


Ruttenetal. (1992) 




i 


75° 


Rutten et al. (1992) 




D 


330 pc 


Rutten et al. (1992) 




Afacc 


~ 10-80 Mq yr-1 


Rutten et al. (1992) 




i^l 


~200 kms-l 


Kaitchuck et al. (1983) 


UXUMa 


AfwD 


0.47 ± 0.07 Mq 


Baptista et al. (1995) 






0.78 ± 0.13 Afg 


Vande Putte et al. (2003) 




M2 


0.47 ± 0.10 Af0 


Baptista et al. (1991) 






0.47 ± 0.07 AfQ 


Vande Putte et al. (2003) 




i 


71.0 ±0?6 


Baptista etal (1991) 




D 


345 ± 34 pc 


Baptista et al. (1991) 




Afacc 


10-8.o±o.2 ^g yr-l 


Baptista et al. (1995) 




Ki 


~200 kms-l 


Schleeel et al. (1983) 



clean continuum in the observations (in contrast to the com- 
plex FUV region). For RW Tri, the continuum fits resulted 
in an accretion rate of 8.3 x 10"^ Mq yr~^ (see Figure O. 
Accounting for the uncertainties in the inclination (75ig °), 

the distance (341^31 pc) and the WD mass (here, an uncer- 
tainty of one-tenth of the solar mass was assumed), the valid 
range for the accretion rate would be 2.7 x 10^^ Mq yr^^ 
to 2.6 X 10~* Mq yr^^, corresponding to a maximum effec- 
tive temperature of the disk between Tcff = 4.9 x 10"' K and 
6.9 X 10-* K. 

The most determined effort to establi sh the s ystem param- 
eters for UX UMa was made by Baptis ta et alJ (IT995) . They 
used a (different) set of high time resolution, low spectral res- 
olution HST spectra, and modeled the time-dependent contri- 
bution of the disk and WD to the continuum to identify the 
phases associated with the eclipse of the WD. From this, they 
inferred a WD mass of 0.47 ± 0.07 A/q, and a mass ratio 
(Afa/MwD) of 1±0.1 and an inclination of 71 ± 1°. They 
then derived a distance of 345 ±34 pc from eclipse maps of 
the system obtained from the combination of UV and i?-band 
light curve s. Since then, most workers, including most re- 
cently Lin nell et al.l (2008), have adopted these parameters 
for UX UMa. There are various reasons to question the ac- 
curacy of the measurements. First, if correct, the WD mass 
in UX UMa is much lower than that of most WDs (includ- 
ing most WDs in CV systems), and the mass ratio is at the 
limit allowed by conserva tive mass transfer. S econd, at the 
temperature determined by 'Baptist a et al.l(ll995h for the WD, 
the WD should have bee n evident in the FUSE spectra of 
UX UMa, but F roning et al. (2003), analyzing time-resolved 
FUSE observations of UX UMa, were unable to identify the 
WD ingress and egress in the light curv e. Despite these con - 
cerns, we too adopt the parameters of Baptista et al.l (Il995h . 
since no alternative values have been proposed with greater 
credibility. They are also listed in Table |2] The value for the 
K v elocity of UX UMa, a gain ^ 200 kms^'^ was adopted 
from lSchlegefetai] (IT981 . 

According to the continuum fit (see Figure|4]i, the accretion 
rate for UX UMa is a little higher than RW Tri. It was deter- 
mined as 1.2 X 10~^ Mq yr"'^ but values within the range of 
6.9 X 10"^ Mq yr-i to 1.9 x 10^^ Mq yr^^ would still be 
consistent with the system parameters considering the uncer- 
tainties in the WD mass, the inclination, and distance. This 
corresponds to a maximum disk effective temperature range 
of Toff = 4.6 X 10* K to 4.9 x 10"* K. 

In order to study the wind geometry and the ionization 



structure of the mass outflow in UX UMa and RW Tri, we 
used PYTHON to develop "point" models for both systems. To 
avoid the computationally demanding brute-force method of 
computing grids of models and searching for the best-match 
spectrum, we adopted a hierarchical approach to searching 
the wind parameters. Our basic methodology was to come 
as close as possible to the observations by adjusting the over- 
all ionization structure in the wind first. For that we started 
with the KWD95 descriptions and the wind parameters de- 
rived bv'Knigge & Drewl (11997) . applied to both our subjects 
of study, RW Tri and UX UMa, and explored different val- 
ues of the mass accretion and the mass-loss rate — these are 
the primary parameters that determine the ionization condi- 
tions. However, we restricted the mass accretion rate to the 
values compatible with the results of disk continuum fits (see 
above). After exploring the A/wind x Afacc parameter space, 
the influence of each of the remaining free wind parameters 
on the spectrum was tested, allowing us to successively obtain 
closer agreement with the observations. When using SV93 
point models, the same approach was used except that the 

Afwind X A/acc spacc was not re-explored and wind param- 
eters which roughly correspond to the best KWD95 model 
were used as a starting point. With this approach we were 
able to explore the sensitivity of the spectrum to the main pa- 
rameters although it did not lead to an exhaustive investiga- 
tion of the complete parameter space. In general, determining 
our "best" models involved subjective criteria since perfect 
matches were never obtained. In these choices, we prioritized 
the overall strength and shape of the major emission lines. 

Applying these strategies, we were able to find good point 
models for both systems using both wind treatments, SV93 
and KWD95. The following sections are dedicated to present- 
ing these models (both their successes and failures) and the 
conclusions drawn from them. Throughout our discussion, we 
focus on the three prominent emission lines of N V, Si IV, and 
C IV. As we shall show, these lines form across a wide range 
of conditions in the wind and, together, probe a significant 
fraction of the outflowing material. In addition, as resonance 
lines in Li- and Na-like ions, they involve relatively simple 
atomic physics which should be well described by our radia- 
tive transfer simulations. We will, however, also comment on 
other spectroscopic features formed by lower ionization state 
metal ions (which exist predominantly in the coolest/densest 
parts of the models) and recombination line features of both 
H and He. 

5. RESULTS 

We now present our point models for both RW Tri and 
UX UMa (see below) and then discuss the similarities and 
differences we infer for the wind structures in these systems. 

Since we allowed for some flexibility in the mass accretion 
rate for a given set of system parameters, we have also al- 
lowed some freedom in the distance we adopt for both objects. 
Specifically, we always renormalized the model spectra to the 
observations in the discussion below. For this renormaliza- 
tion the continua of model and data were compared between 
Si IV and C IV. However, the implied deviations from the dis- 
tances given in Table |2] are well within the uncertainties since 
the point models for both systems have accretion rates very 
close to the values determined by our disk continuum fits. We 
also corrected all pre-eclipse model spectra for the K veloc- 
ity (redshift) of the WD. This shift however was rather small 
— from 0.3 A at the blue end to 0.6 A at the red end of the 
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Table 3 

Parameters for Our Best KWD95 Model for RW Tri and UX UMa. 



MWD 

RWD 
^disk 



. 7 M( 
0.6 A-/^ 
S-0 X : 
2-4 X : 

75°0 
349 pc 
8.0 X : 

5.7 X : 

4.8 X : 

20 
0.1 
3.0 

4.0 X : 

4.5 



0.47 Mq 

0.47 Mq 

9.7 X 10^ cm 
Ti 2.9 X 10^"^ cm 

71 °0 

374 pc 
Afgyr"! 1.3X10-^ A-Jq yr" ^ 2 . 7 X lO^'* Afgyr-l 

5.1 X 10* K 
Mq yr~ll.O X 10~^ Mq yr^^l.O X 10^^ Mq yr^^ 



0.5 
3.0 
4.9 : 
1.0 



15 

0.5 

3.0 

4.9 X 10^ 

4.5 



Notes. 

'^ System and wind parameters from Knigge & Drew! 119971) . These 
parameters were used as a starting point when deriving our point models. 

Distance implied by the renormalization of the model continuum. 
'^ Maximum effective temperature in the disk determined via Equation (T) 

from A/wD, ^WD. AJdisk- 
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Figure 6. Comparison between RW Tri HST observations and our best 
KWD95 model. 



spectrum. 



5.1. RWTri 



Table [3] gives the parameters of our best model for 
RW Tri obtained with the KWD95 wind description via 
our strategy described in Section ID Exploration of the 
-Mwind X Ma_cc spacc for this object led us to prefer Mace — 
8 X 10^^ Mq yr^^ combined with ^fwind — 4.8 x 
10-1° Mq yr-i (- 0.06A/acc)- This best value of Afacc is 
comfortably within the allowed range discussed in Section |4] 
and implies a distance to RW Tri of 13 = 350 pc which is rea- 
sonable. We found that a low value of the mass-loss exponent 
a was required to achieve an acceptable match to the relative 
line strengths in this system. A low value of a means that 
the density near the rotational axis of the system is small, al- 
lowing for a mean higher ionization state in this region. This 
was required in order to produce a sufficiently strong N V 
emission line, relative to C IV. Our exploration of parameters 
indicated that the coUimation of the RW Tri wind is not very 
strongly constrained by the spectra. Only very low collima- 
tions (low values of rf < 10) can clearly be ruled out since 
they lead to clear P-Cygni profiles in C IV which are not ob- 
served. In our best model, a slightly higher collimation was 
adopted than in Knigge & Drew ( 1997), since we found that 
it led to a slightly better match with the data. 

The spectrum obtained for this model at an inclination an- 
gle of 75° is compared to the observations in Figure |6l De- 
tails of the line profiles for the major resonance transitions 
of N V, Si IV, and C IV are shown in Figure |7] These 
are in good agreement with the observations with regard to 
both line width and strength. Since they form under different 
ionization conditions, they are present in quite different re- 
gions of the wind; see Figure[8] In particular, N V probes the 
highest ionization state material which exists in the innermost 
streamlines where the radiation temperature is high owing to 
the abundance of ultraviolet photons radiated by the inner re- 
gions of the disk (see Figure|9l). In contrast. Si IV is confined 
to a relatively narrow equatorial band lying across the base of 
flow streamlines that rise further out in the disk. This region 
has lower Tr since it is exposed to the radiation field from 
the cooler outer disk regions. This, combined with the high 
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Figure 7. Detailed comparison of line profiles for RW Tri and our KWD95 
model during pre-eclipse (left) and mid-eclipse (right). 

electron densities (tIc) at the base of the flow, allows relatively 
low-ionization state material to dominate this region. C IV is 
present across a much wider region of the wind than either 
N V or Si IV, occupying most of the space between the polar 
regions (rich in N V) and the equatorial band of Si IV. C IV 
also exists throughout much of the outer wind allowing for a 
significantly extended region of line formation. This is a con- 
sequence of the relatively uniform distribution of Tr and Ug 
in the outer parts of the flow (see Figure|9l). 

In Figure[8] contours are drawn to indicate the line-forming 
regions. Specifically, these contour lines indicate where the 
last scattering events for Monte Carlo quanta which contribute 
to specific resonance lines took place. As one would expect 
from the ion distribution, C IV has a somewhat extended re- 
gion of line formation, while N V and Si IV are more confined 
to the inner parts of the flow, being concentrated in the polar 
or more equatorial regions, respectively. Given these signif- 
icant differences in region of line formation, it is a powerful 
confirmation of the wind model that the relative strengths and 
shapes of these main line features can be reproduced, both in 
and out of eclipse, by a single model in which the ionization 
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Table 4 

Parameters for Our Best SV93 Models for RW Tri and UX UMa. 




10 20 30 40 10 20 30 40 10 20 30 40 

X [lif cm] 

Figure 8. Number densities for various ions as a function of position in our 
wind model for RW Tri (adopting the KWD95 prescription). Only the inner 
portion of the wind that produces substantial scattering is shown (our com- 
putational grid extended to 10^'' cm). The black contour lines indicate the 
line-forming regions of the main thi'ee resonance lines (upper panels). Specif- 
ically, they are contours of the number of relevant last scattering events in a 
computational grid cell normalized by the cell volume. The four contours are 
drawn for one decade intervals in this quantity which, as expected, decreases 
outward. 



Parameter 


RWTri 


UXUMa 


MwD 


0.7 Mq 


0.47 Mq 


Ma 


0.6 Mq 


0.47 Mq 


i?WD 


8.0 X 10* cm 


9.7 X 10* cm 


^disk 


2.4 X 10^" cm 


2.9 X lO^o cm 


i 


75?0 


71'?0 


£,a 


345 pc 


374 pc 


Mdisk 


8.0 X 10-9 Mq yr-l 


1.3 X 10-* Mq yr-l 


Teff*^ 


5.7 X 10'' K 


5.1 X 10'' K 


A/wind 


4.8 X 10-10 Mq yr-1 


1.0 X 10-9 Mq yr-l 


^min 


RwD 


RwD 


^max 


-Rdisk 


Rdisk 


^min 


2'?9 


1?4 


C'max 


56°3 


36?9 


7 


1 


1 


A 


0.6 


-0.5 


/ 


3 


4 


Rs 


2.0 X 10^1 cm 


4.9 X lOlO cm 


a 


2.2 


1.0 



Notes. 

'^ Distance implied by the renormalization of the model continuum. 

^ Maximum effective temperature in the disk determined via Equation Jl) 

from Af-yvD, i?WD. Mdisij. 
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Figure 9. Temperature, dilution factor, and density structure of the RW Tii 
KWD95 model. 



is also predicted by the model. For He II, the recombination 
line strength is too small in the model but for H I Lya it is 
significantly too strong in the pre-eclipse spectrum. This line 
forms rather deep in the wind models, meaning that its emis- 
sion line flux is significantly affected during eclipse (compare 
model line strengths in Figure |6]l. It is therefore much more 
sensitive to the structure at the base of the wind than the other 
strong emission lines and may indicate a shortcoming of our 
wind models in describing those regions. 

To establish the sensitivity of our conclusions about the line 
formation regions to the choice of how the wind is parameter- 
ized, we have also developed a wind model for RW Tri using 
the SV93 wind prescription. The parameters for this model 
are given in Table |4] and its spectrum is shown in Figure [TO] 
These show comparably good agreement in the spectral fea- 
tures as found with the KWD95 model. Moreover, despite the 
difference in parameterization, they indicate characteristically 
similar geometries for the distribution of ionization states as 
found with the KWD95 wind prescription. 



conditions are controlled by the radiation produced by the ac- 
cretion disk. 

In addition to the three major lines, several metal lines asso- 
ciated with lower ionization state material appear in the spec- 
tra. In the pre-eclipse spectrum, absorption by C II A1333 
and Si III AA1295 — 1303 are present in both model and ob- 
served spectra although Si III is rather too weak in the model. 
C in A1176 also appears in both model and observation but 
its character is less well represented — in particular it appears 
as absorption in the model pre-eclipse spectrum but emission 
in the data. During eclipse this feature is in emission in the 
model but it is too weak. Taken together, this suggests that the 
wind model may place C III too close to the disk plane — a 
better match might be obtained if more of this ion existed in 
the extended (uneclipsed) wind region with less lying along 
our relatively high-inclination line of sight. Recombination 
line emission in both H I Lyman a A1215 and He II A1640 



5.2. UXUMa 

We employed a similar methodology to obtain point models 
for UX UMa with both the KWD95 and S V93 wind prescrip- 
tions. The parameters of these models are given in Tables [3] 
andg] 

Based on the fits to the continuum, we require a lower ef- 
fective temperature for UX UMa meaning that the inner disk 
radiation is redder than in RW Tri. Thanks to this, a good 
match to the relative strengths of emission in N V, C IV and 
Si IV was obtained without having to adopt a low value for a 
(cf. discussion of N V in RW Tri above). An important differ- 
ence in the observed spectra of UX UMa compared to RW Tri, 
however, is the appearance of narrow absorption dips in the 
Si IV profile out of eclipse. We could not reproduce these 
dips via simple changes in either A/^ind or a. To explain 
these features, it is necessary to introduce Si IV ions along 
the observer's line of sight which requires a more equatori- 
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Figure 10. Details of line profiles for the RW Tri SV93 model during pre- 
eclipse (left) and mid-eclipse (right). 

ally concentrated Si IV population. By exploring the velocity- 
law parameters, we found that this could readily be explained 
by a more-rapidly accelerating flow (specifically, we reduce 
the acceleration exponent (3; qualitatively similar results can 
be obtained with a smaller acceleration length Rs). A more 
rapid acceleration reduces the density allowing Si IV to re- 
place Si III around the base of the wind. Together with a rather 
high collimation of the wind, which was necessary to ensure 
the right width and blueshift of the Si IV absorption dips, this 
approach led to a reasonable match of the three main line pro- 
files (see Fi gures [TTI and [TZt . Our low value of /3 = 1.0 is the 
most substantial diff erence in the wind para meters between 
this study and that of'Knigge & Drew! (Il997h : they favored a 
significantly more gradual acceleration (/3 = 4.5) but mod- 
eled only the C IV profile and did not compute the wind ion- 
ization state in detail. Admittedly, they obtained a somewhat 
better fit to the C IV profile than we have (see their Figure 4); 
however, we find that the additional constraints on the wind 
properties by the complex Si IV line profile require a modi- 
fication to the density structure around the base of the wind 
which are consistent with much more rapid acceleration. A 
negative consequence of the change in the density structure 
around the base of the wind is that this model significantly un- 
derpredicts the strengths of features associated with the lowest 
ionization material (e.g., C II and Si III wind features are en- 
tirely absent from the model spectra, while the recombination 
H I Lyman a and He II A1640 lines are too weak). 

While the major emission lines in the model match the data 
well, there is a discrepancy between the slope of the model 
and the observed disk continuum in the pre-eclipse state. The 
model continuum, which reflects that of a steady-state disk, 
is steeper than the observed spectrum, a fact that has been 
known about UX UMa for some time (see, e.g.lK nigge et alJ 
119981) . One possibility is that the run of temperatures in the 
disk simply departs from tha t of a stead y-state disk, an idea 
that was explored in detail bv lLinnell et al. (2008) without, in 
our opinion, finding an obviously better model for the disk 
emission. Our models would suggest that, despite the ap- 
parent mismatch in the continuum slope, the ionization con- 
ditions in the wind required to produce the correct relative 
strengths of the emission lines can be obtained from a stan- 
dard, steady-state disk configuration: if we were to signifi- 
cantly alter the distribution of Teff for the disk in our fiducial 



Figure 11. Comparison between UX UMa HST observations and our best 
KWD95 model. 
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Figure 12. Details of fine profiles for the UX UMa KWD95 mode! during 
pre-eclipse (left) and mid-eclipse (right). 

model, the ionization conditions (and therefore line emission) 
in the wind would change. Although this might be accounted 
for by modifications to other wind parameters, on balance we 
favor an alternative possibility, namely, that the mismatch of 
the continuum slope is the result of occultation of the inner 
portion of the disk along our line of sight. This is certainly 
not an established fact, but UX UMa is known to exhibit pre- 
eclipse dips in the UV that are attributed to this effect; it is 
possible that portions of the disk are occulted all of the time, 
as is known to b e the case in some SW Sex stars, including 
DW UMa (Knigge et al.ll2004 . It seems less likely to us that 
the disk is truncated because this would necessarily reduce 
the peak value of Tcs in the disk and therefore predict a less 
ionized wind. 

The ionization structure of the UX UMa wind model 
(shown in Figure \T3[ is different in detail but similar in char- 
acter to that in the models fit to RW Tri. The main differ- 
ence is the extension of the low-ionization region near the disk 
which is much smaller in this model. This is a consequence 
of the reduction in (3 that was required to reproduce the Si IV 
the absorption dips. 
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Figure 13. Ionization structure of the UX UMa KWD95 model. Again the 
line-forming regions are shown by contour lines (cf. Figure[8j. 
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Figure 14. Details of line profiles for the UX UMa SV93 model during pre- 
echpse (left) and mid-echpse (right). 

As for RW Tri, we were able to transform the wind pa- 
rameters to the SV93 prescription and obtain a wind model 
which produced comparably good spectra to the best KWD95 
model (FigureO. Some minor differences appear in the pro- 
file shapes, in particular the Si IV absorption dips are slightly 
deeper However, once again, the characteristic wind structure 
is qualitatively similar, as are the line-forming regions. 

5.3. Similarities and Differences between the RW Tri and UX 
UMa wind Models 

Compared to previous studies in which the ionization state 
of the wind has been imposed, our approach has the impor- 
tant advantage that the ionization conditions in the wind are 
treated self-consistently with the radiation field. This has di- 
rectly led us to wind models with simple, physically plausible 
ion distributions that are characteristically similar for both the 
objects we have studied with both wind prescriptions. The fa- 
vored picture always involves a polar concentration of N V, 
equatorial Si IV (and lower ions), and C IV filling most of 
the space between. This basic configuration was able to re- 
produce the observations, despite the differences in the UX 
UMa and RW Tri observed spectra. As one would expect, the 



major line formation occurs around the regions where the rel- 
evant ions are most populated and is concentrated within the 
innermost few tens of WD radii. 

Our studies of the two systems have favored mass loss rates 
which are a few percent of the accretion rates (6% and 8%). 
Our models suggest that these values are fairly robust — sim- 
ple changes in the mass-loss rate by a factor of 2 (in either 
direction) leads to results which are strongly inconsistent with 
the observations. 

The main differences between our models for RW Tri and 
UX UMa lie in the details of the wind structure and are con- 
nected to the differences in the KWD95 parameters a (mass 
loss exponent) and j3 (acceleration exponent). In RW Tri, a 
small value of a was necessary to avoid under-producing N V. 
The a-parameter is less strongly constrained for our UX UMa 
point model -a = 0.5 agrees very well with the observations 
and leads to a simple wind structure but higher or lower values 
cannot be strongly excluded. Thus, our models do not rule out 
the possibility that both systems have somewhat similar radial 
dependence of the mass-loss rate per unit area. 

On the other hand, our results relating to the wind accelera- 
tion law do seem to suggest real differences between the sys- 
tems. The /3-parameter cannot be very significantly reduced 
for RW Tri because this would eventually introduce strong 
absorption features in Si IV and C IV, neither of which is ob- 
served. (The high acceleration exponent ensures that the main 
Si IV region is located far above the disk, outside the line of 
sight.) A large value of [3 also allows for a low-ionization 
state region close to the disk (below the main Si IV concen- 
tration) where absorption lines seen in the pre-eclipse spectra 
(e.g., C II) can form. However, the UX UMa spectra also 
place relatively strong constraints on the wind velocity law 
for that system. Here, Si IV absorption dips are seen, favoring 
a higher ionization state closer to the disk compared to RW 
Tri. In our model, this increase in the ionization state was 
realized via a rapid acceleration (low acceleration exponent), 
leading to a low density. We note that th i s is in contrast to 
the study of UX UMa bv lKnigge & Drew! (Il997h in which a 
high /3-value (i.e., slow acceleration) was invoked to account 
for absorption dips in C IV. In their approach, a constant ion- 
ization fraction was assumed meaning that deeper absorption 
arises from higher densities. In our case, we have to account 
for the fact that the density affects not only the optical depth 
directly but also the ionization balance. The net effect of this 
is that lower densities are needed to (see Figures |9] and [TsT i 
obtain the Si IV dips because high densities lead to too low an 
ionization state. Thus, we find that significantly faster accel- 
eration is required to model the UX UMa spectra than those 
of RW Tri, likely pointing to a real physical difference in the 
manner in which the outflows of the two systems are driven. 

A downside of the reduced density around the base of the 
flow in our UX UMa model is that the low-ionization state 
material is pushed very close to the disk. Thus, we lose 
any trace of emission or absorption signatures from very low- 
ionization state material in the model spectra. This suggests 
that a more complex wind structure (e.g., clumping) may be 
needed at the highest inclination to allow for the simultaneous 
appearance of absorption features in both relatively high ion- 
ization state lines (e.g.. Si IV) and significantly less ionized 
material (e.g., C II). We generally found that these classes of 
absorption features could not be simultaneously created with 
the kinematic wind models we have considered. 

An additional, less significant difference between the point 
models of UX UMa and RW Tri is the degree of wind colli- 
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mation. As noted above, there is a certain freedom for this 
parameter in RW Tri since the spectra are not very strongly 
affected by it. We can, however, exclude very low collimation 
since it produces P-Cygni C IV profiles which are in conflict 
with observation. In LfX UMa, winds with very low collima- 
tion can also be excluded for the same reason. But in this 
system, the absorption features in Si IV give additional con- 
straints on the collimation. A moderately high wind collima- 
tion gives the correct blueshift and width of these features, 
introduced by the slowly accelerating wind. 

6. CONCLUSIONS 

We have presented point models for the winds in the eclips- 
ing CVs RW Tri and UX UMa which are able to reproduce 
the main characteristics of the major UV resonance lines of 
N V, Si IV, and C IV both in and out of eclipse. The ion- 
ization conditions in our models are primarily determined by 
the temperature distribution in the wind which, in turn, is de- 
termined by the effective temperature distribution of the disk. 
Accounting for this, we found that agreement with observa- 
tion could be obtained for winds with a characteristic struc- 
ture in which N V forms around the inner edge of the wind 
(exposed to the hottest disk radiation), while Si IV is restricted 
to a region around the base of the flow that extends along the 
equatorial edge of the flow. The equatorial location of Si IV 
allows for narrow absorption dips to appear in the line pro- 
file, as observed in UX UMa out of eclipse. The formation of 
the resonance lines of both N V and Si IV is restricted to the 
innermost boundaries of the wind (within a few tens of WD 
radii of the central object). The outer wind has ionization state 
intermediate between these edge zones and it is here that the 
strong C IV line is formed. The wind volume reasonable for 
this line is significantly more extended than for either of the 
other major resonance lines but still concentrated within ^ 50 
WD radii of the central object. 

Although our focus has been the major resonance lines, our 
models include some wind regions of even lower ionization 
state material which are able to introduce some absorption 
features in the RW Tri pre-eclipse spectra due to, e.g., C II 
and Si III. In agreement with observations, our models also 
predict weak C III emission during eclipse and some recombi- 
nation emission from both H I and He II. The details of these 
features are less well reproduced than those of the main res- 
onance lines; in particular, H I Lya is significantly overpre- 



dicted by our models for RW Tri while He II is dramatically 
underpredicted for UX UMa. The failures are most likely as- 
sociated with the manner in which our wind models describe 
the very base of the wind where it is densest and, perhaps, has 
the most complex flow geometry. 

We have used two different classes of simply parameterized 
flow models (those due to KWD95 and SV93) and confirmed 
that both parameterization are able to reproduce the observa- 
tions with comparable accuracy and that they lead to similar 
conclusions relating to the ionization and structure of the flow. 

The wind structure our models suggest for both RW Tri and 
UX UMa are characteristically similar. Both require wind 
mass-loss rates which are ^6% - 8% of the accretion rate. We 
also found some evidence suggesting that obtaining a wind 
able to simultaneously account for all three of the major reso- 
nance lines may favor models in which the mass-loss per unit 
area of the disk is not strongly coupled to the local disk ef- 
fective temperature. In particular, for RW Tri we found that 
adopting a = 1 for the KWD95 prescription (meaning that 
the mass loss is proportional to the local disk luminosity) sig- 
nificantly underpredicts the N V emission. To reproduce the 
observed N V line strength, we found that a much lower value 
of a ~ 0.1 was required. Qualitatively similar results were 
obtained from the SV93 prescription. In UX UMa, we found 
that the relative line strengths were less constraining on the 
distribution of mass loading in the wind but are also consis- 
tent with a-values somewhat less than 1.0. These conclusions 
are likely to be affected if there are any additional sources of 
ionizing radiation (e.g., a hot boundary layer) that the calcu- 
lations presented here do not include. Nevertheless, they may 
implicate mechanisms other than the radiation force alone in 
launching the wind since it suggests the connection between 
the local mass loss and luminosity in the disk may be weak. 

We found some evidence for real differences in the wind 
properties between the two systems. In particular, our mod- 
els suggest that the acceleration around the base of the flow 
is significantly more rapid in UX UMa than in RW Tri. This 
conclusion mainly stems from the presence of absorption dips 
in the out-of-eclipse Si IV profile of UX UMa which are ab- 
sent from RW Tri. Our attempts to model these features as 
part of the wind argue against a universal kinematic outflow 
model and may suggest that there are real differences in the 
wind launching regions of different systems. However, further 
study of the structure around the base of the wind (including 
consideration of lower ionization state spectral features) is re- 
quired to confirm this and we cannot yet rule out that the ab- 
sorption dips in Si IV might have some different origin (e.g., 
in the rim of the accretion disk). 

Overall, our results support conclusions drawn from previ- 
ous theoretically motivated considerations of the structure of 
CV disk winds a nd their physical origin. In particular, as dis- 
cussed by lDrew & Proga (2000), the mass-loss rates predicted 
from simulations of radiatively driven outflows (Prog a et aP 
119981) are likely too small to account for the strength of ob- 
serve d spectral features. (iDrew & Pro gal 120001; iProga et akl 
120021) . Our results strongly support this conclusion; our wind 
mass-loss rates are indeed significantly higher than those ex- 
pected for line-driven disk wind models with a plausible value 
of the Eddington ratio (see Figure 1 of lDrew & Proga 2000), 
while the disk accretion rates we have used are close to (or 
below) the lower limit of the regime in which radiation pres- 
sure alone might be able to drive a flow. In addition, we have 
concluded that the range of radii from which mass loss is 
launched is significantly more extended than one might ex- 
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pect for a radiative driven flow (i.e., the low a-value which 
we require for RW Tri). This is in line with the discussion of 
iProga et al] ( 1200 2') in which they also argue that line-driven 
disk wind simulations do not succeed in launching fast out- 
flows from as extended a region of the disk as suggested 
by observations. Thus, although the accretion and mass-loss 
rates implied by our models for the eclipsing nova-like vari- 
ables are undoubtedly in the range where radiation forces 
should be dynamically important, additional physics — most 
plausibly magnetic phenomena — Ukely has a part in driving 
the mass loss. 
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